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Abstract 

We calculate the bending angles of light under the strong electric and magnetic fields by a 
charged black hole and a magnetized neutron star according to the nonlinear electrodynamics of 
Euler-Heisenberg interaction. The bending angle of light by the electric field of charged black hole 
is computed from geometric optics and a general formula is derived for light bending valid for any 
orientation of the magnetic dipole. The astronomical significance of the light bending by magnetic 
field of a neutron star is discussed. 
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I. INTRODUCTION 



A nonuniform electric or magnetic background field can induce a continually varying index 
of refraction. At classical level the linearity of the electrodynamics precludes the bending of 
light by electric or magnetic field. Therefore any bending must involve a nonlinear interaction 
from quantum correction. Because of the nonlinear optical properties by Euler-Heisenberg 
interaction jH [2] arising from the box diagram of quantum electrodynamics, the light can 
bend when it passes the neighborhood of an electrically or magnetically charged object. 

Since the maximum available field in a laboratory is of the order B, E/c ~ 10 2 T, it 
is difficult to observe the bending directly in a laboratory [3--7J. An alternative way one 
can think of is observing the effect in the astronomical scale. There are objects, at least 
considered theoretically, with electric or magnetic field strong enough to make such bending 
relevant. One example of light bending by an electrically charged object is when a high 
energy photon passes around a charged black hole with impact parameter greater than the 
Schwarzschild radius [8]. Another example of light bending by a magnetic field is when 
a photon passes the magnetosphere of a magnetized neutron star with extremely strong 
magnetic field of the order 10 8 — 10 11 T. Since neutron stars have a very dense magnetosphere, 
only high energy photons may be observable. Recently the vacuum effect of non-linear 
electrodynamics under strong magnetic field by magnetized stars has been studied widely 

[ME]. 

We will calculate the bending angle of the incident photon under the electric and magnetic 
field. In the previous work of the authors, we set up a simple formalism to find the bending 
angle and the trajectory of light in a Coulombic field at atomic scale [22]. In this paper, we 
consider the bending of a high energy photon when it passes near an astronomical object 
with strong electric or magnetic field. The organization of the paper is as follow. In Sec. II, 
we set up the trajectory equation to calculate the bending angle based on Euler-Heisenberg 
interaction. In Sec. Ill, we consider the bending by the electric field of a spherically symmet- 
ric charged object. In Sec. IV, we consider the bending by the magnetic field of a magnetic 
dipole. We derive a general formula valid for any orientation of the dipole axis and discuss 
some special cases where the dipole axis and the beam direction are aligned in a particular 
way. Finally, in Sec. V, we conclude and discuss the possibility of observation. 
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FIG. 1: Schematic of bending due to nonuniform refractive index. 
II. THE TRAJECTORY EQUATION 

The nonlinear interaction of photons is described by the Euler-Heisenberg Lagrangian 

mm 



4 90m 4 c 
E 2 - c 2 B 2 ) + [(E 2 - c 2 B 2 ) 2 + 7c 2 (E • B) 2 ] . (1) 
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In the presence of the electric field, the correction to the speed of light due to the nonlinear 
interaction is [231433] 

- = 1 - x E) 2 , (2) 

c 45m 4 c 5 V ; ' y ! 

where u denotes the unit vector in the direction of photon propagation, a = 14 for the 
perpendicular mode in which the photon polarization is perpendicular to the plane spanned 
by u and E, and a = 8 for the parallel mode where the polarization is parallel to the plane. 
Throughout the paper all units are in MKS. For magnetic case, E should be replaced by cB 
in the above equation. Because the speed of light depends on the electric field the light ray 
bends in the presence of a nonuniform field. The bending can be calculated by geometric 
optics. The index of refraction due to the background field is given, in the leading order, by 

n=- = l + ^uxE 2 . 3 

v 45m 4 c 5 V ; K J 

The infinitesimal bending of the photon trajectory over Sr can be obtained from the 
Snell's law as 

60 = t&n9— = -|Vn x 5r\ (4) 

n n 



where Sn = Vn • 5v and 9 denotes the angle between u and Vn (see Fig. [I]). Since light 
bends in general toward the direction of greater index of refraction we can write the bending 
in a vector form 

5u = — (<5r x Vn) x u , (5) 
which leads to the trajectory equation 

1 , _ . 

— = — (u x Vn) x u, (6) 

as n 

where s denotes the distance parameter of the light trajectory with ds = \dr\ and 

Since the correction to the index of refraction is generally small in our consideration the 
trajectory equation can be approximated to the leading order as 

= (u x Vn) x u , (8) 

where Uo denotes the initial direction of the incoming photon. Throughout the paper we 
shall assume the photon comes in from x = — oo and moves to +x direction, hence 

uo = (1,0,0), (9) 

and putting Vn = (^1,^2,^3) the trajectory equation becomes 

d 2 x d 2 y d 2 z 

^=°' ^ = ^ 2 ' ds~ 2=m - (10) 
The first equation shows that ds = dx at leading order, which is obviously true, and therefore 
the trajectory equations for y(x) and z(x) in the perpendicular directions to the incoming 
photon are given by 

d 2 y d 2 z 

dx~ 2=r] ^ ^ = % ' (H) 
which will be used in the following analysis. 



III. BENDING BY SPHERICALLY SYMMETRIC CHARGED OBJECT 

We now consider the bending of photon trajectory by a spherically symmetric charged 
object of total charge Q. For this case the bending angle can be calculated in the same way 



4 




FIG. 2: Schematic of light bending in a Coulombic field. 



as in Coulombic case with the electric field 



(12) 
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For a photon trajectory moving to +x axis in the xy plane, choose the unit vector in 
the direction of photon propagation as 

1 



u 



;i,y',o), 



(13) 



V 1 + y' 2 

where prime is the derivative with respect to x. The index of refraction can be written as 

aa 2 h 3 Q 2 (y — xy'^ 2 



n = 1 



(14) 



7207r 2 e m 4 c 5 r 6 (l + y' 2 )' 
For a photon incoming from x = — oo with impact parameter b (see Fig. |2]), the initial 



condition reads 



y(-oo) = 6, Z/'(-oo) = 0, 



and from the first of the trajectory equation (11) 



„ aa 2 Q 2 \\ . y 3y 3 

y =V2 



(15) 



(16) 



3607r 2 eo^c r 6 r 8 ' 

where A e = h/mc is the Compton length of the electron. The total bending angle tp e can be 
obtained by integration 

POO 

y'(oo) = / r] 2 dx = tan ip e ~ ip e (17) 
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aa 2 Q 2 f\ e . 



By putting y = 6 in 772, for the leading order solution, we obtain 

6407re ftc \ 6 

The bending always occurs toward the center of the charged object as in the bending by 
gravitational field. 

To compare the bending by electric field with the bending by gravitation, let us consider 
a charged non-rotating black hole with mass Ai and charge Q. The bending by gravitational 
field is well-known 

Note that, from ip g oc 1/b and ip e oc 1/& 4 , the bending by electric field can be important at 
short distance. The charge and angular momentum per unit mass(J/Ai) is constrained by 
the mass of the black hole, in Planck units, as 



Q 2 + {J/Mf < M 2 . (20) 

For non-rotating (J = 0) charged black hole, restoring the physical constants, the total 
electric charge is constrained by the condition 

47re 

We can parameterize the charge as 



< GM 2 . (21) 



Q = y/4Tre GMZ, (22) 

with < £ < 1. Then the magnitude of the bending angle by electric field can be written as 

_ aa 2 i 2 GM 2 ( \ e \ 4 _ aa 2 ? bcM ( A e \ 4 
LPe ~ 160 he V b ) ~ 640 h \b J Vg ' ^ ' 

To estimate the size of the bending, let us compare the two bending angles for the 
maximally charged (extremal, £ = 1) stellar black hole of ten solar mass M. = 10.M sun = 
2 x 10 31 kg. Since our formalism is based on flat space time, not on general relativity, the 
impact parameter should be large enough. We consider the case when the impact parameter 
is ten times the Schwarzschild radius, b = 10r s h ~ 300 km, at which 

(p g = 1.98 x 10~ x rad; 

<p e = 5.47 x 10 V g = 1.08 x 10 1 rad (for the perpendicular mode a = 14). (24) 



The bending by electric charge dominates the gravitational bending. Even for non-extremal 
charged black hole with £ = 0.1, the electrical bending, (p e = 1.08 x 10 _1 rad, is comparable 
to the gravitational bending. 

This example shows that the bending by electric field can be large. However, it may not 
survive the screening by electron-positron pair creations of the strong electric field. As well 
known, as the electric field strength approaches the QED critical field 

E c = m 2 c 3 /eh = 1.32 x 10 18 V/m (25) 

the electric field is susceptible to electron-positron pair creation. In the region where the 
electric field is of the order or higher than E c , the vacuum is unstable and the electric field is 
highly screened by the pair creation [21 133]. Only photons entering the region with electric 
field below E c can have a chance to be observed. 



To simplify the discussion we may write (23) as 




^«°.°83(£h/M^)(T). (26) 



•10 

where % * s the ratio of the black hole mass to the solar mass, 



Eh is the Coulomb field at radius b: 



x = ^-- (27) 



with Q given by (22), and b c denotes the radius at which E h = E c given by 




& c « 3654/lRo, (29) 



where R is the Schwarzschild radius of the black hole R = 2GM.jc . Eq. (26) shows that 
at b, for which the electric field strength is below E c , the ratio is bounded by 

K < O083 , a x (3Q) 



fx uo. 

which shows that the screening by pair creation severely bounds the valid range of the 
bending angle by a black hole with mass larger than the solar mass. The light bending 
can be comparable to the gravitational bending only for black holes with mass substantially 
smaller than the solar mass. 
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FIG. 3: Schematic of light bending by magnetic dipole. The magnetic dipole is located at the 
origin, the dotted line is the photon trajectory, the dashed line is the projection of outgoing photon 
path on the xz-plane, and 4>h (4> v ) is the bending angle of horizontal (vertical) direction. 

IV. BENDING BY MAGNETIC DIPOLE 

We consider the bending of photon trajectory by a magnetic dipole. Obviously, the bend- 
ing by magnetic dipole should depend on the orientation of dipole relative to the direction 
of the incoming photon. We consider the bending by a magnetic dipole located at the origin 
with arbitrary orientation (see Fig. [3]). We take the direction of the incoming photon as x 
axis, horizontal direction as y axis, and vertical direction as z axis. For the magnetic dipole 
M located at origin, we define the directional cosines a M ■ x. I M ■ y. ~ .\/ • .: such 
that M = MM. The magnetic field by the dipole is given by 




(31) 



and the index of refraction due to this background magnetic field is 




(32) 



Taking the unit vector in the direction of photon propagation as 




1 



(33) 



u = 



the index of refraction can be written as 



n = 1 + 



aa 2 h 3 e ( fi M\ 2 



1 



[{B z y' - B y z'f + {B z - B x z'f + {B y - B x y'f] . 

(34) 



45m 4 c 3 \ 47r J 



1 + y' 2 + z 
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To the leading order, the index of refraction is written explicitly as 



n 



1 + 



aa 2 h 3 e 



45m 4 c 3 \ An J r 
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{my 2 



X 



z 2 ) + 



3axy + 3^yz j 



+ <^ 7(2^ 2 - x 2 - y 2 ) + 3ax2 + 3f3yz 



(35) 



and from the photon trajectory equation (J 1 1 p we have 

2 



„ 2aa 2 ^ 3 e / (J>oM\ ' 
^ ~ 45m 4 c 3 V 4?r J 



r 



1% 



| (/3(2y 2 - a; 2 - ^ 2 ) + 3azy + 3-yyz) (A/3y + 3ax + 3-fz) 

+ (~f{2z 2 - x 2 - y 2 ) + 3axz + 3(3yz) (-2 7 y + 3pz) } 

3axy + 371/2;) 2 



(/3(2y 2 - x 2 - , 2 ) 



+ (7(2 



2 — x 



y 2 ) + 3axz + 3{3yz)' 



(36) 



2aa 2 h 3 e ( HqMX 
A5m 4 c 3 \ An J 



{my 



2 2 2\ 

x — z ) 



lOz 

r 12 



+ 3axy + 3 r yyz) (—2(3z + 3jy) 

+ (^{2z 2 - x 2 - y 2 ) + 3axz + 3(3yz) (4 7 z + 3ax + 3/3?/) j 
j (/3(2y 2 - a; 2 - z 2 ) + 3aa;?/ + 372/2) 2 

+ ( 7 (2z 2 - x 2 - y 2 ) + 3axz + 3(5yz) 2 } 



(37) 



The total bending angle can be obtained by integration by putting y = b and z = in (36) 



and (37), which gives for horizontal ((fh = y'(oo)) and vertical (cp v = z'(oo)) deflections, 

2 i4 



_^^(^ (15q2 + 41 , 2 + 1672) 

5n aa 2 e c ( pL M\ \\ 

3-2 6 h \~^r) ¥ fl ' 



(38) 
(39) 



An important comment is in order. The result is valid as long as the polarization of the 
photon remains, throughout the trajectory, perpendicular or parallel to the plane spanned 
by u and B. It can be easily checked that this condition is met only when the magnetic 
dipole axis is either in z direction or in the xy plane. When the dipole axis is in none of 
these directions the photon polarization remains neither in pure perpendicular nor in pure 
parallel mode, even if it started in one of the modes. In this case the birefringence effect 




FIG. 4: Schematic of light bending by magnetic dipole when the photon path is perpendicular to 
the dipole moment and traveling on the equator of the dipole. 

takes place, and the perpendicular and the parallel modes of the light ray split. Since the 
splitting occurs continually on every branch-out rays, this cascade of splitting results in the 



initial single light ray branching out to a light bundle. The bending angles (38) and (39), 
with a either 14 for the perpendicular mode or 8 for the parallel mode, then provide an 
envelope for the maximal and minimal bending angles of the light bundle in the horizontal 
and perpendicular direction, respectively. 

Let us now consider some special cases which will allow us to compare our result with 
the previous work. First, we consider the case when the photon path is perpendicular to 
the dipole moment and traveling on the equator of the dipole. Assume that the magnetic 
moment directs along z and the incident photon is coming from x = — oo (see Fig. [4]). This 
is the specific case considered by Denisov et al. [llj. Taking a = (5 = and 7 = 1, the 
magnetic field on the dipole equator (xy plane) is given by 



47T r 3 



By symmetry, there is no vertical (z) bending, so ip v = 0, and the horizontal (y) bending 
angle is given by 

7T aa 2 e c ( /i M\ 2 \\ 

^ h = -24^r[-^r) ¥■ (41) 

This result agrees with Denisov et al. (see the Eqs. (4) and (5) in |llj). 

Second, consider the case when the photon path is parallel or anti-parallel to the dipole 
axis. Assume that the dipole at the origin directs along the —x axis in the xy plane (see 
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FIG. 5: Schematic of light bending by a magnetic dipole, located at origin directing — x direction, 
when the photon path is parallel to the dipole 




FIG. 6: Schematic of light bending by magnetic dipole when the photon path is perpendicular to 
dipole moment and passing the axis of the dipole. 

Fig. |5j). Taking a = — 1 and — 7 = 0, the magnetic field on the xy plane is given by 

B = ^± E ((y*-2x 2 )x-3xyy). (42) 

Also there is no bending in z direction, ip v = 0, and the bending angle in y direction is given 
by 

Finally, consider the case when the photon path is perpendicular to dipole moment and 
passes the north or south pole. Locate the dipole at the origin directing along the y axis in 
the xy plane (see Fig. [6]). Taking a = 7 = and (3 = 1, the magnetic field on the xy plane 
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5?r aa 2 e c ( /i M \ 2 



is given by 



B = ^i(3*yx+(2y 2 -a; 2 )y). (44) 
There is also no vertical bending, (p v = 0, and the horizontal bending angle is given by 

4l7r aa 2 e c ( a M\ 2 A 4 

^--3^—{^r) ¥■ (45) 

This configuration gives the maximum possible bending; The gradient of index of refraction 
is maximal along this direction. 

Let us now compare the bending by magnetic field with the gravitational bending. We will 



consider the possible maximum bending given by Eq. (45 ) for a strongly magnetized neutron 
star with solar mass Ai = -M S un = 2 x 10 30 kg and radius ro = 10 km. For extremely strongly 
magnetized neutron stars the magnetic field at the surface can be strong as B s = 10 8 — 10 11 T. 
Parameterizing the impact parameter in units of the radius b = (r with £ > 1, the bending 
by a magnetic field can be written as 

4l7r aa 2 e c -Xi , 

= s^-S-^ (46) 

where we have used B s = /ioM/47rrQ, the magnetic field strength at the neutron star surface. 
For the magnetic field strength of B s = 10 9 T of a neutron star, the maximum value of the 
bending angle is given by 

<p m = 1.40 x HT 4 rad, (47) 

which is for a ray glancing the north or south pole where the bending is maximal. This is 
much smaller than the gravitational bending cp g = 0.59 rad. 

By increasing the surface magnetic field a larger bending angle can be obtained. However, 
as in the bending by charged black hole, the bending is constrained by the screening of 
electron-positron pair creation of strong magnetic field above the critical strength 

B c = = 4.4 x 10 9 T . (48) 

en 



To see this we may write (46) as 



ip g V10/ \RoJ \B t 

where B^ is the magnetic field at radius b given by 



B h = B s (^y, (50) 
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and Rq denotes the Schwarzschild radius Ro = 2GA4/c 2 . To avoid the screening we require 
B b < B c , which gives the bound 



where b c denotes the radius at which = B c . Thus, a larger bending can be obtained with 
large b c /R , but it is obviously bounded by the surface magnetic field strength through the 
relation 



For a magnetar with B s = 100-B C , for example, we get b c = 4.6ro, and considering that 
tq/Rq ~ the magnetic bending is expected to be, at most, a few percent of the 

gravitational bending. 

V. DISCUSSION 

We have studied how photons can be bent when they travel through the strong electric 
or magnetic field of compact object like a charged black hole or a neutron star. We cal- 
culated the bending angles according to the nonlinear electrodynamics of Euler-Heisenberg 
interaction. Our calculation shows that the bending by electrically charged astronomical 
object can be comparable or larger than the gravitational bending. However, the screen- 
ing by the electron-positron pair creation strongly bounds the electric bending so that the 
electric bending can be significant compared with the gravitational bending only for black 
holes with mass substantially smaller than the solar mass. We also found a general formula 
for light bending by magnetic field, valid for any orientation of the magnetic dipole. Our 
calculation shows that for a magnetar with surface magnetic field B s = 10 11 T the magnetic 
bending can be, at most, a few percent of the gravitational bending. 

Since the magnetic bending is expected to be small compared to gravitational bending any 
chance of observation may be realized only when the experiment detects the small variation 
over the gravitational bending. One way to observe the light bending by magnetic field may 
be using the birefringence [9] that the bending of perpendicular polarization is 1.75(=14/8) 
times larger than the bending of parallel polarization. Even in the region where the bending 
by magnetic field is weak compared with the gravitational bending, by eliminating the overall 




(51) 




(52) 
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FIG. 7: Schematic of detecting the light bending by a neutron star and its binary system of equal 
mass. The bending angles at t = and t = T/2 can be different by the contribution of magnetic 
field of the neutron star. 

gravitational bending, the polarization dependence of the bending by magnetic field may be 
tested if the allowed precision is sufficient enough. 

Another possibility is detecting the time variation in bending. The magnetic axis of 
a magnetar is likely to be different from the rotational axis. In this case the magnetic 
bending by a rapidly spinning magnetar will add a small wiggling effect over the gravitational 
bending. 

The neutron star in a binary system [19] may also be used to detect the magnetic bending. 
Most of the neutron stars are isolated stars, less than one hundred are known to be in binary 
systems with nondegenerate stars. Assume for simplicity that the nondegenerate companion 
star has the same mass as the neutron star. If we consider only the bending by gravitation, 
the bending at t — and the bending at a half orbital period later t = T/2 will be the 
same. However, when the bending by magnetic field is included, the bending angles will be 
different by the relative position of the two stars (see Fig. [7]). 
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